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Abstract 

Vascular endothelial growth factor-A (VEGF) is produced by most cancer cells as multiple Isoforms, which display distinct 
biological activities. VEGF plays an undisputed role in tumour growth, vascularisation and metastasis; nevertheless the 
functions of individual isoforms in these processes remain poorly understood. We investigated the effects of three main 
murine isoforms (VEGF188, 164 and 120) on tumour cell behaviour, using a panel of fibrosarcoma cells we developed that 
express them individually under endogenous promoter control. Fibrosarcomas expressing only VEGF188 (fsl88) or wild type 
controls (fswt) were typically mesenchymal, formed ruffles and displayed strong matrix-binding activity. VEGF164- and 
VEGF120-producing cells (fs164 and fsl20 respectively) were less typically mesenchymal, lacked ruffles but formed 
abundant cell-cell contacts. On 3D collagen, fs188 cells remained mesenchymal while fsl64 and fsl20 cells adopted 
rounded/amoeboid and a mix of rounded and elongated morphologies respectively. Consistent with their mesenchymal 
characteristics, fsl88 cells migrated significantly faster than fsl64 or fs120 cells on 2D surfaces while contractility inhibitors 
accelerated fs164 and fsl20 cell migration. VEGF164/VEGF120 expression correlated with faster proliferation rates and lower 
levels of spontaneous apoptosis than VEGF188 expression. Nevertheless, VEGF188 was associated with constitutively active/ 
phosphorylated AKT, ERKl/2 and StatS proteins. Differences in proliferation rates and apoptosis could be explained by 
defective signalling downstream of pAKT to FOXO and GSK3 in fsl88 and fswt cells, which also correlated with p27/p21 
cyclin-dependent kinase inhibitor over-expression. All cells expressed tyrosine kinase VEGF receptors, but these were not 
active/activatable suggesting that inherent differences between the cell lines are governed by endogenous VEGF isoform 
expression through complex interactions that are independent of tyrosine kinase receptor activation. VEGF isoforms are 
emerging as potential biomarkers for anti-VEGF therapies. Our results reveal novel roles of individual isoforms associated 
with cancer growth and metastasis and highlight the importance of understanding their diverse actions. 
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Introduction 

Vascular endothelial growth factor-A (VEGF) plays a funda- 
mental role in tumour growth, vascularisation and metastasis and 
exists as multiple isoforms derived by alternative splicing of the 
VEGF gene [1]. Mouse and human proteins of 120/121, 164/165 
and 188/189 amino acids respectively, represent major VEGF 
splice variants with distinct properties and expression patterns. 
These isoforms dilTer in terms of binding affinities to the 
extracellular matrix and receptor activation. Tumours display 
highly variable levels of relative isoform expression, with VEGF- 
164/165 and VEGF 120/ 121 generally being the most predom- 
inant and VEGF- 188/ 189 relatively less abundant [2]. VEGF 
signals through tyrosine kinase receptors VEGFRl/flt-1, 
VEGFR2/flk-l and VEGF3/flt-4 [3]. VEGF also binds neuropilin 



co-receptors (NRP-1 and NRP-2), which lack tyrosine kinase 
activity but regulate the function of VEGF receptors as well as 
other receptor tyrosine kinases (RTKs) [3] . The dilferent affinities 
to matrix, displayed by the various VEGF splice variants generate 
gradients in vivo and result in different signalling responses, which 
are important for angiogenesis [4,5]. VEGF also has complex 
functions in angiogenesis-independent aspects of tumour growth 
and tumour cells have been shown to express functional VEGF 
receptors [6,7,8] but the role of individual VEGF isoforms in these 
processes remains poorly understood. 

VEGF and its receptors are now major targets of several cancer 
therapies. Anti-VEGF agents such as the humanised neutralising 
anti-VEGF antibody bevacizumab as well as several VEGF 
receptor kinase inhibitors are being used to treat many types of 
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cancer. However, not all patients respond to anti-VEGF therapy 
and therefore biomarkers that can predict clinical response are 
being actively pursued [9]. Indeed, several recent retrospective 
clinical studies have identified the short soluble isoforms of VEGF 
(VEGF 120 and 110) as promising potential biomarkers for 
predicting responsiveness to bevacizumab [10,11,12]. Pancreatic, 
breast and gastric cancer patients with higher circulating levels of 
short VEGF121 and VEGFllO isoforms were shown to have a 
shorter median overall survival but were more likely to benefit 
from treatment with bevacizumab. Importantly, the short VEGF 
isoforms were identified out of a comprehensi\ e range of potential 
biomarker candidates that were examined in these studies [8,9,10]. 
These studies highlight the importance of understanding of the 
complexities associated with the functions of individual VEGF 
isoforms in cancer. 

Immortalization and transformation of mouse embryonic 
fibroblasts results in cell lines with tumorigenic potential. Such 
cells derived from genetically modified animals are useful for 
studying the role of specific genes in cancer progression [13]. 
Using this approach we developed VEGF isoform-specific 
fibrosarcomas [14] by H-rasV12 transformation of embrv'onic 
fibroblasts derived from transgenic animals expressing single 
VEGF isoforms and wild type controls [15,16]. These cells 
produced tumours in vivo that displayed distinct vascular patterns, 
similar to those seen in corresponding transgenic animals during 
development. Furthermore, the tumours displayed differences in 
response to treatment with vascular targeting agents highlighting 
the importance of VEGF isoform expression in treatment outcome 
[14,17]. When grown in vitro, the fibrosarcoma cells exhibited 
major differences in growth rates, with VEGF164 and VEGF120 
expressing cells (fsl64 and fsl20 respectively) proliferating 
significanfly faster than VEGF188 (fsl88) and wild type control 
cells (fswt) [14] suggesting that VEGF isoform expression also 
controls tumour cell growth characteristics. Initial growth rates 
were also faster for fsl64 and fsl20 tumours in vivo, which are 
consistent with the highly proliferative phenotype of the cells [14]. 

Proliferation and survival of cancer cells is maintained by 
constitutive activation of multiple complex signalling pathways 
that ordinarily maintain tight control in normal cells [18]. 
Mutated proto-oncogenes such as Ras can cause aberrant 
signalling through die Raf/MEK/ERKl/2 and PI3K/AKT 
signalling cascades that regulate tumour cell growth and survival. 
Multiple signalling mechanisms also converge to co-ordinate 
tumour cell invasiveness and migration. In particular, assembly of 
the actin cytoskeleton, under the control of the Rho family of 
GTPases, is required so that protrusions and contractile forces can 
be generated to propagate cell movement [19]. Cancer cells move 
either individually or collectively as multicellular clusters and 
maintain the capacity to switch from one mode of movement to 
another. Individual cell migration is classified as either "mesen- 
chymal" or "amoeboid" and each mode has different require- 
ments for contractility, attachment to matrix and proteolysis [20]. 

VEGF is known to contribute to intracrine and autocrine 
tumour cell growth, survival, migration and invasion [7,8]. We 
hypothesise that individual isoforms of VEGF might influence 
various aspects of tumour cell behaviour in a distinct manner. 
Using the fibrosarcoma cells we developed, which express single 
VEGF isoforms and, in parallel, wild type control fibrosarcoma 
cells capable of expressing all isoforms, we investigate here the 
consequences of endogenous isoform expression on tumour cell 
morphology, proliferation, survival and migration. Using this 
system, in which VEGF expression is maintained under the 
control of its endogenous physiological promoter, we have 
discovered that individual isoforms are associated with distinct 



roles in terms of tumour cell behaviour and gained insight into 
signalling pathways that are linked with these effects. 

Materials and Methods 

Ethics statement 

Animal experiments were conducted according to United 
Kingdom Animals (Scientific Procedures) Act 1986 (UK Home 
Office Project Licence PPL40/31 10) and with local University of 
Sheffield ethical approval. 

Reagents 

Blebbistatin, U0126, Y27632 and pyridone 6 (P6) were 
purchased from Calbiochem. SU 11248 (sunitinib) was purchased 
from EC laboratories. Antibodies to integrin fil, ILK, pAKT 
(ser473), p70S6 (thr389), pGSK-3 (ser21/9), pFOXO-1 (ser256), 
p-Stat3 (tyr705), VEGFR2, pVEGFR2 (tyr996), pVEGFR2 
(tyrll75) and GAPDH were from Cell Signalling. Anti- actin, 
pERK, tERK and P-tubuIin antibodies were from Sigma. N- 
cadherin, p21 and p27 antibodies were from BD Biosciences and 
LIE antibody was from R&D Systems. VEGFRl (ab32152 and 
ab2350) and pVEGFRl {lyrl333) antibodies were from Abeam. 
Recombinant VEGF164, VEGF165 and VEGF120 proteins were 
purchased from Peprotech. Recombinant VEGF 188 protein was 
purchased from Reliatech GmbH. 

Cell culture 

Fibrosarcoma cell lines (fsl88, fsl64, fsl20 and fswt) were 
developed as we described before [14], from fibroblasts derived 
from transgenic mouse embryos expressing single VEGF isoforms 
VEGF188, VEGF164, VEGF120 and wild type controls respec- 
tively [15,16]. Cells were maintained in Dulbecco's modified 
Eagle's medium (DMEM), \Q% FCS, 600 ng/ml G418 and 2 |ig/ 
ml puromycin (Invitrogen, UK). The mouse endothelial H5V cell 
Une was a gift from Dr Annunciata Vecchi [2 1] . H5V cells were 
maintained in DMEM containing 10% FCS. Human umbilical 
vein endothelial cells (HUVEC) from pooled donors were 
purchased from PromoCell and were maintained in endothelial 
cell growth medium (PromoCell). 

Preparation of conditioned media 

Cells were grown to 50% confluence and then switched to 
serum-free DMEM for 48 h. Conditioned media were concen- 
trated 30-fold with Microcon concentrators (Merck MiUipore) and 
normalised to cell numbers before they were analysed by 
immunoblotting. 

Preparation of collagen, fibronectin and laminin matrices 

A thin coating of monomeric collagen was prepared by 
incubating culture plates for 2 h with 50 |J,g/ml rat tail collagen 
type. I (BD Biosciences) diluted in 0.02 N acetic acid. Bovine 
fibronectin (Invitrogen, UK) or laminin- 1 (Sigma) were diluted to 
5 |Xg/ml in PBS and used to coat dishes for 2 h. To prepare a thick 
layer of fibrillar matrix, ice-cold rat-taU type I collagen was diluted 
to 1.5 mg/ml with lOXDMEM and sterile dH20. The solution 
was neutralised with NaOH and allowed to polymerise at 37°C 
before plating the cells on top. 

Adhesion and cell spreading assays 

Adhesion assays were performed in 96-well plates pre-coated 
with fibronectin, laminin or collagen and blocked in 0. 1 % BSA. 
CeUs (3x10*) suspended in serum-free DMEM/0.1% BSA were 
plated in each well, and plates were immediately incubated at 
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37°C for 45 min. Adherent cells were fixed, stained with 0.5% 
crystal violet and lysed in 2% SDS. Absorbance was measured at 
570 nm using a BMG FLUORStar Galaxy microplate reader. For 

monitoring spreading, cells were re-suspended in DMEM/0.1% 
BSA and plated in Lab-Tek 4-well permanox chamber slides 
(Nalge-Nunc) (1x10' per well) pre-coated with collagen, laminin 
or fibronectin. The slides were incubated at 37°C for 30 min to 
3 h and then cells were stained as described below. 

Analysis of cell morphology and cytoskeletal F-actIn 
staining 

Live cells were stained with 2.5 |J.g/ml CellMask plasma 
membrane stain (Invitrogen, UK) for 5 min, fixed in 'i.T'/o 
formalin and mounted in Vectashield (Vector Laboratories). 
Analysis of F-actin was performed as we described previously 
[22]. Fluorescence images were taken with a Leica DMI4000B 
fluorescence microscope and using LASAF control and analysis 
software. 

Growth curves and doubling times 

Cells were plated in 6-well plates (4x10* per well). Viable cells 
were counted daily for 4-6 days using a Vi-Cell automated cell 
viability analyzer (Beckman Coulter). Media were replenished 
every other day. Doubling times were established by plotting a 
graph of log 2 (N/NO) (NO = initial cell numbers and N = end- 
point average cell numbers) against time followed by calculating 
the inverse slope of the hnear part of the curve. 

Colony formation in soft agar 

Cells were suspended in medium containing 0.33% agarose and 
plated over a 0.5% agarose underlay in 6-well plates (2x10^ cells/ 
well). Three weeks later, colonies were stained with 0.005% crystal 
violet and counted in 10 random fields using a 4x objective. 

Analysis of apoptosis 

Apoptotic cell death was measured by flow cytometry using the 
ApoAlert MitoSensor Kit (Clontec) which detects alterations in 
mitochondrial membrane potential. Apoptosis was also quantified 
using the Cell Death Detection ELISA'^'"'^^ kit (Roche Diagno- 
sitics), which measures cytoplasmic DNA-histone nucleosome 
complexes. For both assays cells were plated at a density of 10* 
cells/cm^ and apoptosis was determined 48 h later. Apoptosis in 
tumour sections was detected by TUNEL assay using the ApoTag 
Plus Peroxidase In Situ Apoptosis Detection kit (Chemicon 
International). Apoptotic cells were scored using a 40 x objective 
in 6-9 random fields of view per tumour section avoiding areas of 
necrosis. Data are expressed as mean number of TUNEL positive 
apoptotic cells per mm^. 

Migration assay 

Migration was measured using Ibidi cell culture inserts. Cells 
(3.5x10*) in 70 |xl medium were plated into each insert 
compartment and grown for 48 h. The inserts were removed to 
create a gap/wound of 500 |j,m into which the cells migrated. The 
cells were imaged immediately and at intervals for up to 24 h with 
a Nikon Eclipse phase contrast microscope equipped with a Digital 
Slight DS camera and NIS-Elements software. Wound closure was 
quantified with an automatic image analysis algorithm [23] using 
the CAIMAN image analysis website (http://www.caiman.org. 
uk/). 



Subcutaneous tumour generation 

Fibrosarcoma cells (1x10^ in 50 (jl serum-free DMEM) were 
injected subcutaneously into the rear dorsum of female 8-12 week- 
old SCID mice. Tumours were excised when they reached a mean 
diameter of 6-8 mm (usually 10-14 days post-implantation) and 
were either embedded in paraffin and sectioned for analysis of 
apoptosis, or were frozen for protein extraction. 

Immunoblotting 

Proteins from cells and tumours were extra('t(;d with NP40 Lysis 
Bufier (Biosource, Invitrogen) supplemented with phosphatasi; and 
protease inhibitors. Equal amounts of protein (10-50 |ig/lane) 
were separated on NuPAGE Novex gels (Invitrogen), transferred 
to nitrocellulose or PVDF membranes and immunoreactive bands 
were visualized using ECL reagents (GE Healthcare). 

Statistical analysis 

Data were analysed by using a one-way analysis of variance 
(ANOVA), followed by the Tukey-Kramer post-test test for 
multiple comparisons using GraphPad Prism software for Mac 
OS X. A 2-way ANOVA followed by a Bonferroni post-test was 
used for analyzing multiple groups with more than one variable 
(treatment and cell type). In all cases, difierences between groups 
were described as significant if the probability was <0.05. 

Results 

Fibrosarcoma cells expressing individual VEGF isoforms 
differ in morphology and adhesion properties 

Fsl88 as well as fswt cells, which express all isoforms including 
VEGF188, display typical mesenchymal ft-aturcs (Figure la,b). 
The cells are elongated and spindli'-shaj)ed, have extended 
processes and ruffles and do not form (■xtcnsi\(- cell-cell 
associations. In contrast, fsl64 do not display a typical mesenchy- 
mal morphology. These cells grow in very close association with 
each other and closely align longitudinally forming long multicel- 
lular chains. Fsl20 cells also grow in close association with each 
other, although less so than fsl64 cells, and also have fewer 
extended processes than fsl88 and fswt cells. On fibrillar collagen, 
fswt and fsl88 cells retained their mesenchymal features while 
fsl64 cells were rounded/ amoeboid and fsl20 cells displayed a 
mix of rounded and elongated features (Figure 2a). All fibrosar- 
coma cells stained diffusely for F-actin (Figure 2b), a characteristic 
feature of Ras transformation, which elevates ERK MAPK 
signalling and uncouples Rho-GTP signalling from stress fibre 
formation [24,25]. Inhibition of ERK activator MEK with U0126 
reversed the effects of transformation on F-actin, and re- 
established stress fibers in all the cell fines (Figure 2c). In addition, 
fsl20 cells assumed a more mesenchymal morphology and 
displayed more extended processes. Fsl64 cells expressed stress 
fibres upon MEK inhibition although their distinct close ceU-ceU 
aUgnment was still evident, suggesting that this feature was not 
MEK-dependent. 

Fswt and fsl88 cells adhered better to uncoated plastic and 
collagen than fsl64 and fsl20 cells (Figure 3a). On the other hand, 
fsl64 and fsl20 cells adhered better to laminin than fswt and fsl88 
cells. All four cell lines adhered to and also spread avidly and 
rapidly on fibronectin (Figure 3a, b). All cells also spread on 
laminin but fsi64 and fsl20 cells failed to spread on collagen, even 
after 3 h and in the presence of serum (Figure 3b). AU the 
fibrosarcoma cells expressed integrin (31 and integrin-linked kinase 
(ILK) (Figure 3c), an adaptor protein that is recruited to integrin 
fSl cytoplasmic domains and associates with actin at focal 
adhesions [26]. ILK expression levels were higher in fswt and 
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fswt fs188 fs164 fs120 



a) 



b) 



Figure 1 . Morphology of fibrosarcoma celis expressing single VEGF isoforms. a): Phase contrast images of cells grown on uncoated tissue 
culture plastic, b): Cells grown on plastic and stained live with CellMask Orange. Scale bars, 50 (im. 
doi:1 0.1 371/journal.pone.01 0401 S.gOOl 



fsl88 compared to fsl64 and fsl20 cells. All fibrosarcomas also 
expressed N-cadherin (Figure 3c), a classic mesenchymal marker 
[27]. Levels of N-cadherin were up-regulated in fsl64 and fsl20 
cells, which also formed more cell-cell contacts (also see Figure 1). 



Fibrosarcoma cells expressing individual VEGF isoforms 
proliferate at different rates and display differences in 

levels of spontaneous apoptosis 

We previously showed that fsl64 and fsl20 cells proliferated on 
plastic at significandy faster rates than fswt and fsl88 cells [14]. 
Integrin-matrix interactions play an important role in controlling 
cell proliferation. However, neither collagen nor fibronectin had 
any significant effect on growth, and the differences in prolifer- 



fswt 



fs188 



fs164 



fs120 




Figure 2. iVIorphology of fibrosarcoma cells on thick fibrillar collagen and after inhibition of MEK. a): Cells on a thick layer of type I 
collagen; b): cells grown on uncoated plastic and stained for F-actin. c): Cells treated with IVIEK inhibitor U0126 (10 (ilVI) for 18 h and stained for F- 
actln. Scale bars, 50 nm. 
doi:1 0.1 371/journal.pone.01 0401 5.g002 
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Figure 3. Fibrosarcoma cell adiiesion and spreading, a): Cell adiiesion to uncoated plastic or plastic coated with collagen I, fibronectin or 
laminin-1 at 45 min. * and ** represent significant differences, (*p<0.05 and **p<0.01) by two way ANOVA followed by Bonferroni post-test, b): Cell 
spreading on fibronectin and laminin at 30 min, on a thin layer of collagen at 3 h or on collagen for 3 h in the presence of 2% FCS (collagen+S). Cells 
were stained live with with CellMask orange. Scale bars, 50 jim. c): Immunoblot analysis for pi integrin (P-ltg) ILK and N-cadherin. Blots were 
normalised to actin. Results are representative of at least 3 independent experiments. 
dor:1 0.1 371 /journal.pone.01 0401 5.g003 



ation between the cell lines persisted (Table I). Growth in soft 
agar, an indication of tumourigenic capacity, was also significantly 
different between the cell Unes; fsl64 and fsl20 cells producing 
significantly more colonies than fswt and fsl88 cells (Table II and 
also see Figure SI). Apoptosis was significantly more pronounced 
in fswt and fsl88 cells than in fsl64 and fsl20 cells (Figure 4a,b). 
Similarly, in vivo, there were significandy more apoptotic cells 
within viable regions of fsl88 solid tumours compared to tumours 
generated by fsl64 and fsl20 cells (Figure 4c,d). 

Fibrosarcoma cells expressing individual VEGF isoforms 
migrate at different rates on 2D surfaces 

Fswt and fsl88 cells migrated significantly faster than fsl64 and 
fsl20 cells on 2D uncoated plastic (Figure 5). The faster migration 
of VEGF188-expressing cells was consistent with their mesenchy- 
mal morphology and presence of ruHles (see Figure 2b), an 
indication of intense migratory activity [28]. Migration rates on 
thin collagen were similar to migration rates on uncoated plastic 
(data not shown). 



Migration of fsl64 and fs120 cells is actinomyosin 
contractility-dependent 

Cell morpholog}' and migration are dependent on the 
cytoskeleton and actinomyosin contractility, to provide traction 
force in mesenchymal movement, and cortical contraction for 
rounded amoeboid movement [19,29]. Cells were treated with 
pharmacological inhibitors of contractility, including ROCK 
inhibitor Y27632, myosin II ATPase inhibitor blebbistatin and 
JAK inhibitor P6. The latter was shown to inhibit actinomyosin 
contractility in melanoma cells [30]. Figure 6a shows that all three 
inhibitors caused fsl64 and fsl20 cells to switch to a mesenchymal 
morphology (compare Figure 6a with Figure 1) and accelerated 
their migration but not that of fswt and fsl88 cells (Figure 6b). P6 
inhibited MLC phosph()r\iation in all the fibrosarcoma cells 
(Figure 6c) as prc\iously shown for melanoma cells [30] but 
interestingly also induced a marked up-regulation of integrin pi 
(Figure 6c). 
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Table 1. Population doubling times (hours) of fibrosarcoma cells grown on 2D surfaces. 









fswt 


fs188 


fs164 


fs120 


Plastic 


18.36±0.29 (n = 6) 


17.89±0.336 (n = 6) 


***13.89±0.57 (n = 6) 


•**13.67±0.26 (n = 6) 


Collagen 


16.98±0.92 (n = 6) 


17.17±0.24 (n = 3) 


••13.15±0.09 (n = 3) 


*14.49±0.46 (n = 3) 


Fibronectin 


19.41±1.175 (n = 6) 


18.87±1.55 (n = 3) 


"•13.55±0.24 (n = 3) 


***14.59±0.69 (n = 3) 



Population doubling times (hours) of cells grown on plastic, collagen or fibronectin are means of 3-6 independent experiments ± SEM. 

*p<0.05, 

**p<0.01, 

***p<0.001 values represent differences between fswt/fs188 cells versus fs164/fs120 cells {two way ANOVA followed by Bonferroni post-test). 
doi:l 0.1 371 /journal.pone.Ol 0401 S.tOOl 



Fibrosarconna cells expressing individual VEGF isoforms 
exhibit differences in growth and survival-associated 
signalling 

We next investigated signalling mechanisms that could poten- 
tially account for the observed difTerences in cell growth, survival 
and migration. Signalling via AKT promotes growth through 
activation of the mTOR 1 complex and downstream targets such 
as p70S6 kinase involved in cell growth and Gl cell cycle 
progression [18]. AKT and p70S6 were more highly phosphor- 
ylated in fswt and fsl88 cells, compared to fsl64 and fsl20 cells 
(Figure 7a) and this was independent of cell density (data no 
shown). Overnight exposure to PI3K inhibitor LY29004 reduced 
constitutive phospho-AKT and p70S6 levels (Figure 7a). AKT 
also signals through FOXO and glycogen synthase kinase-3 (GSK- 
3) to regulate survival, metabolism and growth [31]. GSK-3 is 
active when unstimulated and becomes switched off upon 
phosphorylation by AKT. FOXOl is also phosphorylated and 
inactivated by AKT and this action blocks its transcriptional 
activity on genes that promote apoptosis and regulate metabolism. 
Paradoxically, levels of pGSK-3 and pFOXOl did not parallel 
AKT and p70S6 phosphorylation (Figure 7b). Higher levels of 
phosphorylated FOXOl were detected in fsl64 and fsl20 cells in 
which basal pAKT was low, suggesting that AKT signalling was 
uncoupled from this downstream branch of the pathway. 
Similarly, pGSK-3 was not upregulated as might be expected in 
VEGF188-expressing ceUs. Both FOXOl and GSK-3 mediate ceU 
cycle arrest by regulating expression and stability of cyclin- 
dependent kinase inhibitors (CDKIs) [32,33]. In fsl88 cells, CDKI 
p27 was significantly upregulated, while in fswt cells expression of 
CDKI p21 was higher than in the other lines (Figure 7b). 
Increased phosphorylation of p44 and p42 ERKs was also seen 
in fswt//fsl88 versus fsl64/fsl20 cells (Figure 7c). Constitutive 
activation ofJAK/Stat3 is a feature of cancer cells where its major 
function is to promote proliferation and survival [34]. Phospho- 
Stat3 levels were consistently higher in fswt and fsl88 cells, 
compared to fsl64 and fsl20 cells; pan-JAK inhibitor P6 blocked 
Stat3 phosphorylation in all cell lines (Figure 7d). A similar pattern 
of pStat3 was evident in sub-confluent as in post-confluent cultures 
(data not shown), which demonstrated that unlike other systems its 



expression was not density-dependent [35]. Stat3 phosphorylation 
was also up-regulated in fsl88 tumours in vivo (Figure 7e). 
Although Stat3 signalling is primarily associated with proliferation 
and survival, previous studies have shown that sustained activation 
of Ras/Raf/ERK led to the production of leukemia inhibitory 
factor (LIF), which then signaled through JAK/Stat3 in an 
autocrine or paracrine manner to mediate cell cycle arrest [36]. 
Figure 7f shows that fsl64 and fsl20 cultures produced signifi- 
cantly more secreted LIF than fswt and fsl88 cells. It is therefore 
unlikely that signalling through LIF was responsible for the 
differences in levels of active JAK/Stat3 that were seen between 
the cell lines. 

Fibrosarcoma cells express tyrosine kinase VEGF 
receptors, which cannot be activated by VEGF 

We next explored the possibility that the observed differences in 
growth and survival-associated signalling were due to VEGF 
isoform-dependent autocrine actions. All fibrosarcoma cells 
expressed similar levels of VEGF receptor 2 (VEGFR2/flk-l) 
(Fignre 8a) which could potentially mediate autocrine signalling. 
However, VEGFR2 was not constitutively phosphorylated/acti- 
vated in any of the cells (not shown), and furthermore, exogenously 
added recombinant murine VEGF 164 failed to activate its 
phosphorylation either at tyrosine 996 (Figure 8a) or at tyrosine 
1175 (Figure 8b). Downstream ERKl/2, AKT and Stat3 were 
also not further activated/phosphorylated beyond basal levels by 
recombinant VEGF, results that were consistent with the failure of 
VEGF to activate its receptor(s) (Figure 8c). In parallel, VEGF 
induced robust VEGFR2 and ERKl/2 phosphorylation in 
HUVEC and H5V mouse endothelial cells (Figure 8a,b). In 
addition to VEGF 164, other recombinant VEGF isoforms 
(VEGF120, VEGF188 and VEGF165) were tested over a range 
of concentrations (1, 10 and 100 ng/ml) but in all cases induction 
of VEGFR2 and/or ERKl/2 phosphorylation was not detected 
(data not shown). The VEGF receptor tyrosine kinase inhibitor 
SU 11248 [37] blocked both basal and VEGF-induced VEGFR2 
and pERK phosphorylation in H5V cells (Figure 8a) and HUVEC 
(data not shown) but interestingly had no effect on basal levels of 
pERKl/2 in the fibrosarcomas (Figure 8d). These results suggest 



Table 2. Colonies formed by fibrosarcoma cells grown in soft agar. 









fswt 


fs188 


fsie4 


fs120 


Soft Agar 


23.4±2.85 (n = 3) 


24.98±4.4 (n = 3) 


*33.8±3.7 (n = 3) 


*35.8±4.2 (n = 3) 



Number of colonies in agar are means of three independent experiments ± SEIVl. 

*p<0.05 represents significantly more colonies formed by fs164 or fs120cells compared to fswt or fs188 cells (ANOVA followed by Tukey-Kramer post-test). 
doi:l 0.1 371 /journal.pone.Ol 0401 5.t002 
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Figure 4. Spontaneous apoptosis in fibrosarcoma cells and solid tumours, a): Nucleosome generation in cells measured by ELISA at 48 h. 
OD values are averages of 5 independent experiments ± SEM. b): Changes in cell mitochondrial potential at 48 h. Values are averages ± SEM, from 4 
independent experiments. *p<0.05, **p<0.01 represent significance using repeated measures ANOVA followed by Tukey-Kramer post-test, c): 
Apoptosis in tumour sections determined by TUNEL assay. Arrows indicate typical apoptotic cell nuclei. Bars, 50 |im. d): IVIean number of TUNEL 
positive apoptotic cells ± SEM per mm^ tumour section (viable regions only); n = 5-10 tumours per group. *p<0.05, **p<0.01 and ****p<0.0001 
represent significance using repeated measures ANOVA followed by Tukey-Kramer post-test. 
doi:1 0.1 371 /journal.pone.01 0401 5.g004 
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Figure 5. Migration of fibrosarcoma cells expressing single VEGF isoforms. a): Wound closure was measured at intervals up to 24 h. Each 
point represents data obtained from 2-5 independent experiments ± SEM; b): representative culture images of the four fibrosarcoma cell lines 
showing wound closure at 18 h. 
doi:1 0.1 371 /journal.pone.01 0401 5.g005 
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Figure 6. Analysis of morphology and migration of fibrosarcoma cells exposed to contractility inhibitors, a): Cells were treated with P6 
(5 |a,M), blebbistatin (5 \iM) or Y27632 (10 \iM) for 18 h. Representative images of P6 treated cells by phase contrast microscopy; blebbistatin-treated 
cells stained for F-actin; Y27632 treated cells stained with CellMask Orange. Scale bars, 50 ^m. b): Migration was quantified by measuring wound 
closure. Images were collected at 12 h for the fswt and fsISS cells, at 18 h for the fs120 cells and at 24 h for the fs164 cells to correlate with 40-60% 
wound closure in untreated conditions. Contractility inhibitors were added at the same concentrations used in (a) 18 h before the start of the assay 
and replaced in fresh media at assay start. Each point represents data obtained from 3 independent experiments ± SEM, each conducted with 2-4 
replicates. *p<0.05, ***p<0.001 values represent significance over untreated cultures (ANOVA followed by Tukey-Kramer post-test), c): Cells were 
incubated overnight with P6 (5 ^IVI). Equal amounts of proteins (30 ng) were analysed for integrin (31 (|3-ltg) and pIVlLC. A representative immunoblot 
is shown. 

doi:1 0.1 371/journal.pone.01 0401 5.g006 



that although the fibrosarcoma cells express VEGFR2, this 
nevertheless remains inactive/inactivatable. Furthermore, these 
data strongly indicate that the markedly higher levels of 
constitutively active pERKl/2, which were evident in fswt/fsl88 
compared to fsl64/fsl20 cells were unlikely to have arisen 
through VEGF-dependent autocrine activity through tyrosine 
kinase VEGF receptors. However, SU11248 caused a significant 
reduction in constitutive levels of AKT and StatS phosphorylation 
of in aU the fibrosarcoma cells (Figure Bd); after overnight 
incubation with SU 11248, the phosphorylation levels of both 
AKT and StatS returned back to their original basal levels (data 
not shown). 

The fibrosarcoma cells also expressed full length VEGFRl/fltl 
as well as a truncated variant most likely consisting of the 
intracellular domain of the receptor (see supplementary data. 
Figure S2) previously described in endothelial and breast cancer 
cells [38]. The truncated variant was constitutively phosphorylated 
at tyrl333 in the fibrosarcomas cells but this phosphorylation 
could not be downregulated by SU 11248. Phosphorylated 
VEGFRl was also evident in H5V ceUs; in these ceUs SU11248 



caused a small reduction in VEGFRl phosphorylation (see Figure 
S2). 

Discussion 

In this study, we demonstrate that endogenous single VEGF 
isoform expression is linked to distinct growth, survival and 
migratory characteristics in Ras transformed fibrosarcoma cells. 
Our model system allowed us to investigate three individual 
isoforms of murine VEGF (VEGF188, VEGF164 and VEGF120), 
when expressed under the control of the gene's physiological 
promoter, and hence has a distinct advantage compared to over- 
expression systems that may alter normal endogenous signalling 
interactions. We show that VEGF188 expression is linked to a 
slower rate of tumour cell proliferation and decreased survival. 
Furthermore, VEGF188 expression was associated with strong 
cell-matrix matrix interactions and a mesenchymal morphology 
and mode of cell motility in vitro. In contrast, VEGF164 and 
VEGF 120 endogenous expression was associated with more rapid 
cell proliferation, increased sur\dval and a rounded amoeboid 
morphology and mode of motility in vitro (summarised in 
Figure 9). Dilferences between the various cell phenotypes are 
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Figure 7. Signalling pathway analysis in fibrosarcomas expressing single VEGF isoforms. a-d): Where indicated cells were incubated 
overnight with JAK inhibitor P6 (5 |iM), or PI3K inhibitor LY20990 (10 |iM). Equal amounts of proteins (10-30 |ig) were analysed for a); phospho-AKT 
(p-AKT) and p70S6; b): pGSKB, pFOXOI, p27 and p21; c): phospho-ERKs1/2 (p-ERK) and total ERK (t-ERK); d): phospho-Stat3 (p-Stat3). e): Proteins 
extracted from solid tumours (50 |ag/lane) and analysed for pStat3. All blots were normalized with actin or GAPDH and are representative of at least 
three independent experiments, f): Concentrated conditioned media normalised against cell numbers, analysed for LIF expression. 
doi:1 0.1 371/journal.pone.01 0401 5.g007 



very unlikely to have arisen by clonal selection during the 
generation of the cell lines, as several individual clones isolated 
after transformation maintained identical properties to their 
corresponding parental lines (data not shown). 

The cells we developed are of mesenchymal origin and 
accordingly they all constitutively expressed classic mesenchymal 
markers such as N-cadherin (Figure 3c), alpha smooth muscle 
actin and vimentin (unpubhshed data). However, despite not 
displaying typical mesenchymal features, in terms of their 
morphology and migratory characteristics, the fsl64 and fsl20 
cells express higher levels of N-cadherin compared to the more 
typically mesenchymal fswt and fsl88 cells. We postulate that the 
higher levels of N-cadherin protein facilitate and sustain the 
distinct cell-cell contacts that fsl64 and fsl20 cells form in 2D 
culture. Indeed, previously it was shown that upregulated N- 
cadherin expression mediated cell-cell adhesion by localising to 
cell-cell contacts in invasive epithelial MDCK cells after under- 
going epithelial-to-mesenchymal transition (EMT) [39]. 

DiflFerences in adhesion to various matrix proteins are also 
suggestive of dilferences in integrin expression and/or status of 
activation. These differences are likely to involve both alpha and 
beta integrin subunits. While aU fibrosarcomas expressed integrin 
P 1 , fs 1 88 cells expressed higher levels of ILK. ILK stabilises (3 1 
adhesions, is associated with integrin engagement/ activation and 
was found to be over-expressed in mesenchymal cancer cells 
[40,41]. Therefore, upregulated ILK expression could potentially 
contribute to a more adhesive mesenchymal fsl88 phenotype. 
Amoeboid morphology, which is less dependent on integrin- 
matrix interactions, is more readily seen in culture when cells are 



grown on fibrillar collagen matrices. Indeed, fsl64 and to a large 
extent fsl20 cells were rounded on collagen (Figure 2a). Interest- 
ingly, fsl64 and fsl20 cells have rounded cell morphologies when 
grown as subcutaneous tumours in vivo and fsl88 tumours retain 
an elongated mesenchymal cellular morphology (see Figure 4c). 

The extent of anchorage independent growth in agar was 
consistent with growth on solid surfaces. Growth in agar reflects 
tumorigenic potential and results in vitro correlate with our 
previous findings in vivo, where expression of VEGF 164 and 
VEGF120 was associated with a more rapid initiation of tumour 
growth [14]. 

Although VEGF is a pro-survival factor, there are suggestions 
that VEGF 188 may elicit pro-apoptotic signals as shown in 
Figure 4. Chondrocytes from homozygous VEGF188/188 mice 
apoptosed at high rates during development, and could be rescued 
by exogenous VEGF 164 [42]. In addition, over-expression of 
VEGF 189 in breast cancer cells induced apoptosis via NRP-1 
[43]. Interestingly, the fswt and fsl88 cells markedly over-express 
NRPl compared to the other two cell lines (see Figure S2d). It is 
therefore possible that NRP-1 may regulate pro-apoptotic 
processes associated with VEGF 188 expression in the fibrosarco- 
mas. Our results are therefore in agreement with the above studies 
and suggest that increased pro-apoptotic signalling is an intrinsic 
property of VEGF188 expression. It was, however, intriguing that 
the pro-survival and pro-proliferative signalling proteins AKT, 
ERKl/2 and Stat3, were more highly activated in VEGF 188- 
expressing cells. A potential explanation for this discrepancy came 
from the analysis of downstream PI3K/AKT effectors. While 
AKT target p70S6 was phosphorylated, unexpectedly, FOXOl, 
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Figure 8. Analysis of VEGFR2 expression/activation in fibrosarcoma cells. a,b,c): Fibrosarcoma cells, HUVEC or H5V endothelial cells were 
stimulated with 100 ng/ml recombinant VEGF164 for 10 min. Equal amounts of proteins (10-30 |jg) were analysed for a): Phospho-VEGFR2 (tyr996), 
total VEGFR2 and pERK1/2. SU1 1248 (10 ^M) was added 5 min prior to VEGF; b): phospho-VEGFR2 (tyr1175) and total VEGFR2± 10 ^iM SU1 1248; c): p- 
ERK1/2, pAKT and pStat3. d): Cells were incubated with SU1 1248 (10 |iM) for 1 h before proteins were extracted and analysed for pERK1/2, pAKTand 
p-Stat3 (p-Stat3). All blots were normalized with an antibody to p-tubulin and are representative of at least three independent experiments. 
doi:1 0.1 371/journal.pone.01 0401 5.g008 



another key target of AKT and regulator of cell cycle and pro- 
survival signalling [44] , was more highly phosphorylated in fs 1 64 
and fsl20 cells (Figure 7b). Inactive/phosphorylated FOXOl 
correlated with lower levels of p21 and p27 in fsl64 and fsl20 
cells. In contrast, in fsl88 cells, p27 was markedly upregulated and 
in wild type cells, p2 1 was more highly expressed, thus correlating 
with their slower proliferation rates. The mechanism(s) through 
which active AKT signalling couples to p70S6 but not FOXO and 
GSK-3 in fsl88 remains unclear. PP2A phosphatase regulates the 
activity of FOXOl and GSK-3 [45,46]. It is possible that high 
phosphatase activity dephosphorylates FOXO and GSK3 in fsl88 
cells. Several studies have suggested that the pro-survival function 
of AKT can be overridden [47] . Chronic AKT activation can 
contribute to apoptosis induction, while many stress-inducing 
agents and chemotherapeutic drugs activate AKT. Indeed, pro- 
apoptotic agents often activate survival signals, including Stat3 and 
ERKs in addition to AKT, as part of a protective response [48] . 
The fact that Stat3 was highly active in both fsl88 cells and solid 
tumours in vivo suggest that high levels of Stat3 activity represent 
an intrinsic property consequent to VEGF 188 expression. 

High levels of constitutively phosphorylated AKT and presence 
of lameUipodia suggests a link between the PI3K/ AKT pathway 
and the mesenchymal mode of migration of fswt and fs 1 88 cells. 
Our results are in agreement with studies showing that mesen- 
chymal cell lines exhibit high levels of AKT phosphorylation 
[41,49]. PI3K/AKT signalling is also a critical component of the 
epithelial-mesenchymal transition [50,51]. Ras-mediated activa- 
tion of ERKs promotes activation of Rac and lameUipodia 
formation to drive tumour cell motility [52]. Our results are 



therefore in line with these studies and suggest that the features of 
cells expressing VEGF- 188 are linked with activated Ras/PI3K 
and Ras/ERK. 

Differential requirement for integrin engagement and matrix 
adhesion as well as actinomyosin contractility are major features 
that distinguish amoeboid from mesenchymal migration 
[20,53,54] . In mesenchymal cells, Rho/ROCK driven contractil- 
ity supports the formation of focal adhesions and contractile stress 
fibers that can generate traction force for forward movement, 
while Rac and Cd42 are required for the extension of protrusions 
at the cell's leading front [28,55]. Mesenchymal movement is also 
dependent on integrin-matrix interactions. On the other hand, in 
rounded amoeboid cells integrin-matrix adhesive interactions are 
weak, while cortical myosin II generated contraction is essential for 
migration in 3D matrices [53]. Several studies have shown that 
loss of contractility switches cells over to a mesenchymal mode of 
movement [30,54,55]. Inhibitors of Rho/ROCK or contractility 
impaired amoeboid movement but had minimal effects on 
migration in mesenchymal tumour cells thus illustrating that 
Rho/ROCK signalling/ contractility are essential in amoeboid but 
dispensable in mesenchymal motility [30,54,55] . When exposed to 
contractility inhibitors, the fsl64 and fsl20 cells acquired 
mesenchymal characteristics (they spread more and formed more 
protrusions) and were able to migrate faster on 2D surfaces. Our 
results are in agreement with the studies described above in that 
inhibitors of contractility had no significant effect on fswt or fs 1 88 
cells but altered the morphology and motility of fsl64 and fsl20 
cells. In our study, all the fibrosarcomas expressed similar levels of 
phosphorylated MLC, a read-out of actinomyosin contractility 
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Figure 9. Schematic of proposed signalling interactions and morphological and motility characteristics of fibrosarcomas 
expressing single VEGF isoforms. Differences in PI3K/AKT, Ras/Raf/MEK/ERK and JAK/Stat3 signalling between fswt/fs188 and fs164/fs120 cells 
are shown. Upregulation/downregulation of activities or levels of expression of individual proteins are indicated by up or down arrows. VEGF188 
expression is associated with reduced proliferation rates, increased apoptosis and typical integrin-dependent mesenchymal morphology and 
migration mode. VEGF164/VEGF120 expression is associated with faster proliferation rates, increased survival and rounded/amoeboid integrin- 
independent morphology and mode of motility. Migration of fs164/fsl20 but not fswt/fsl88 cells is contractility-dependent. 
doi:1 0.1 371 /journal.pone.01 0401 5.g009 



(Figure 6c), suggesting tliat contractility alone is not sufficient to 
dictate morphology and migration. Inhibition of Jak/Stat3 
signalling abrogated the phosphorylation of MLC as described 
previously for melanoma cells [30] thus establishing this pathway 
as a major driver of contractility also in fibrosarcomas. However, 
while in melanoma cells, pStatS correlated with amoeboid mode of 
movement, in our system pStatS expression correlated with a 
mesenchymal morphology that was independent of contractility 
for migration. Others have shown that mesenchymal hepatocel- 
lular cancer cells expressed higher levels of pStatS compared to 
their epithelial counterparts [41]. Interestingly, inhibition of JAK 
signalling resulted in a marked upregulation of integrin fil 
expression, which correlated with a switch to a more mesenchymal 
morphology, particularly in fsl64 and fsl20 cells. 

We hypothesised that the observed differences in the signaUing 
characteristics of the fibrosarcoma cells could have arisen through 
isoform-dependent autocrine activity, mediated by tyrosine kinase 
VEGF receptors. However, although the fibrosarcoma cells 
expressed RTK VEGF receptors, these could not be activated/ 
phosphorylated by recombinant VEGF isoforms. Exogenously 
added VEGF also failed to activate downstream signalling to 



pERKl/2, pAKT and pStatS in the fibrosarcoma cells, in marked 
contrast to its robust activation of VEGF signalling in endothelial 
cells. Furthermore, adding back different recombinant VEGF 
isoforms had no effect on fibrosarcoma cell proliferation (supple- 
mentary Figure S3) and had no influence on cell morphology (data 
not shown). Sunitinib (SU11248) a potent RTK inhibitor with 
VEGFR selectivity [37] had no effect on basal pERKl/2 levels 
further confirming that the intrinsic pERKl/2 levels were 
VEGFR-independent. Sunitinib did inhibit pAKT and pStatS 
but this is most likely due to inhibition of other RTK targets of 
sunitinib, such as PDGFRfi. Previously it was shown that VEGF 
produced by renal carcinoma cells promoted tumour cell growth 
through activating NRP-1 receptors and downstream Ras [7]. 
While our results cannot exclude an involvement of NRP-1 in 
autocrine signalling, the fact that various recombinant VEGF 
isoforms failed to activate pERK, pAKT or pStat3 signalling or 
alter the proliferation of the fibrosarcoma cells does not support 
such an involvement. 

Taken together our data suggest that the inherent differences in 
signalling, proliferation, survival and migration between our 
different cell lines are a result of complex interactions governed 
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by endogenous VEGF isoform expression, that are independent of 
simple ligand-receptor activation. VEGF 189 has been shown to 
localise to the nucleus [56,57], which could be significant here. 
Although the function of VEGF in the nucleus has not been 
established, nuclear localised VEGFD was found associated with c- 
myc and RNA polymerase II and regulated c-myc-dependent gene 
transcription in fibroblasts in a VEGF receptor-independent 
manner [58]. It is therefore possible that VEGFA isoforms also 
have intracrine activity that is not dependent on classic receptor 
activation. 

In summary, we present, for the first time, evidence that 
VEGF 188 expressed by tumour cells is associated with mesen- 
chymal motility, slower proliferation, and increased apoptosis, 
while VEGF 164 and VEGF 120 are associated with a more 
rounded/amoeboid morphology and increased proliferation and 
survival. Both amoeboid and mesenchymal modes of movement 
are important in metastasis and different targeting approaches 
may be needed to attack different invasive strategies. In this 
respect, it is of interest that VEGF 189 isoform over-expression was 
shown to correlate with metastasis and poor prognosis in colon 
and lung cancer [59,60]. Deciphering the various processes by 
which individual VEGF isoforms contribute to tumour growth is 
important for the design of cancer therapies. Until now, such 
therapies focused on targeting VEGF signalling with an aim to 
impair angiogenesis. The demonstration that other non-angiogen- 
esis dependent aspects of tumour growth are reliant on VEGF 
isoform expression adds to the complexity of VEGF signalling 
within the tumour microenvironment and necessitates a better 
understanding of the molecular mechanisms involved. 

Supporting Information 

Figure SI Growth of fibrosarcoma cells in soft agar. 

Colonies formed by fibrosarcoma cells grown in soft agar and 

imaged using a lOx objective. 

(TIF) 

Figure S2 VEGFRl and NRP-1 expression by fibrosar- 
coma cells, a) Schematic diagram showing VEGFRl domains 
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